Crossing the phantom divide without phantom matter 
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A class of braneworld models can lead to phantom-like acceleration of the late universe, but 
without the need for any phantom matter. In the simplest models, the universe contains only cold 
dark matter and a cosmological constant. We generalize these models by introducing a quintessence 
field. The new feature in our models is that quintessence leads to a crossing of the phantom divide, 
w — — 1. This is a purely gravitational effect, and there is no phantom instability. Furthermore, the 
Hubble parameter is always decreasing, and there is no big rip singularity in the future. 



I. INTRODUCTION 

Observations of supernovae redshifts, cosmic mi- 
crowave background anisotropies and the large-scale 
structure provide increasingly strong evidence that the 
late-time universe is accelerating. If general relativity is 
a correct description of large-scale gravity, then the ac- 
celeration typically originates from a dark energy field 
with w = p/p < — i. The simplest model is LCDM, 
where the dark energy is the vacuum energy (w = —1). 
It gives a very good fit to the data Q, but with an un- 
naturally small and fine-tuned value of A. Quintessence 
scalar fields, with w > — 1, produce more general dynami- 
cal behaviour, but they do not improve the fit to the data, 
and also do not solve the theoretical problems faced by 
LCDM. "Phantom" scalar fields, with w < — 1, have the 
same theoretical problems as quintessence, but in addi- 
tion, classically they have the unphysical behaviour that 
the energy density increases with the expansion of the 
universe, and they also lead to instability of the quan- 
tum vacuum. 

Current data does not rule out a phantom value, w < 
— 1. The WMAP 3-year data, in combination with large- 
scale structure and SN data, allows for w < — 1 ; in a flat 
general relativistic model Q, 



w = -1.04 ± .06. 



(1) 



Models with w < — 1 but without the problem of quan- 
tum instability are therefore still an interesting possibil- 
ity. An effective w such that w < — 1 can occur in mod- 
ified gravity theories, without any phantom matter that 
renders the quantum vacuum unstable. Instead, grav- 
ity itself produces phantom-like acceleration. Sahni and 
Shtanov 3] showed that a class of braneworld models ex- 
hibited effective phantom behaviour. In these models, 
the 4-dimensional brane universe contains only matter 
and a cosmological constant A. A 5-dimensional gravita- 
tional effect screens A, leading to phantom dynamics. 

These models are a variant of the Dvali-Gabadadze- 
Porrati (DGP) braneworld model, generalized to cosmol- 
ogy by Deffayet 4\. The standard DGP model is a self- 
accelerating model without any form of dark energy, and 



the effective w is always non-phantom. However, there 
is another branch of DGP models, with a different em- 
bedding of the 4D brane universe in the 5D bulk space- 
time. The self-accelerating DGP model is the (+) branch. 
The DGP(— ) model is very different. It does not self- 
accelerate, but requires dark energy on the brane. It ex- 
periences 5D gravitational modifications to its dynamics, 
which effectively screen the dark energy. At late times, 
as gravity leaks off the 4D brane, the dynamics deviates 
from general relativity. The transition from 4D to 5D 
behaviour is governed by a crossover scale r c , the same 
as in the DGP(+) branch. 

The simplest DGP(— ) model has a cosmological con- 
stant A, and we follow Lue and Starkman [f| and call 
this the LDGP model. The dynamics of these models 
were investigated by Lue and Starkman p| , and observa- 
tional constraints on the models have been considered by 
Lazkoz et al. The energy conservation equation for 
CDM remains the same as in general relativity, but the 
Friedman equation is modified: 



p + 3H P = , 



- = -(p + A) 

I c O 



(2) 
(3) 



[The DGP(±) branches have TH/r c in Eq. ©.] These 
equations imply 



H = -- 



1 - 



1 



v/l + 4r2(p + A)/3 



(4) 



Equation J2J shows that at early times, i.e., for H ^> r^ 1 , 
the general relativistic Friedman equation is recovered, 
but at late times, the H/r c term is important and the 
Friedman equation is nonstandard. 

At late times, gravity leakage screens the cosmological 
constant, leading to an effective dark energy 0,0] 



A 1 H 
PeS = A - 3— , 



(5) 



where p e s and u> e ff = p e fi/peS are effective quantities in 
a general relativistic interpretation of LDGP expansion 
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history, i.e., they describe the equivalent general relativ- 
ity model: 



H 2 = ^(p + Pcs), 

p cff + 3H(1 + w cS )p cS = . 



It follows that 



1 + ffldf 



H 



r c Hp cS ' 



(6) 
(7) 

(8) 



and since H < by Eq. @, we have effective phantom 
behaviour, ui c ff < —1, provided that p e ff > 0, i.e., for 
H < r c A/3. 

In general relativity with phantom matter, w c g < — 1 
implies that H eventually becomes positive, i.e., the uni- 
verse eventually super-accelerates, which can lead to a 
"big rip" singularity. By contrast, in LDGP H is al- 
ways negative, and there is no big rip singularity. In 
LDGP however, the phantom divide w c ff = — 1 cannot 
be crossed: w c s is always less than —1 for p c s > 0. At 
higher redshifts, when H > r c A/3, we have w c s > —1, 
but w e g never passes through —1, so there is no crossing 
of w c ff = — 1. The effective phantom picture only holds 
for p eff > 0. 

In this paper we generalize the LDGP model by intro- 
ducing a quintessence field with w q =const (> —1), some- 
times called "quiessence" . We will call this the QDGP 
model. We will show that crossing of the phantom divide 
does occur in QDGP models, unlike the LDGP limit. 
The interesting fact is that in general relativity with a 
dark energy component, phantom crossing cannot occur 
with a single field 7] (unless non-minimal coupling be- 
tween dark matter and dark energy is allowed Q)- In 
the QDGP model, crossing is possible without resorting 
to a complicated model with multiple fields or additional 
non-gravitational interactions. 



II. THE QDGP MODEL 

The Friedman equation is Eq. ©, with 

PcS = Pq - 3— . (9) 

Since w q is constant, the conservation equations give 

p = p Q (l + z) 3 , p q = p q0 (l + z)^ 1+ ^K (10) 
We define dimensionless density parameters, 



fin 



Po 
3H* 



PqO 

3#o 2 



(11) 



Then Eqs. © and © give 



H v 



(12) 




FIG. 1: The deceleration parameter q = —(1 + E/HoE 2 ) 
against redshift, for (Q m ,fl q ,Wq) = (0.28,0.8,-1) [solid 
(black) curve], (0.28,0.8,-0.99) [dotted (red) curve], and 
(0.28,0.8,-0.95) [dashed (blue) curve]. 



so that 



2y/n rc = n m + n q - 1 > o. 



(13) 



Thus the flat QDGP model mimics a closed GR 
quiessence model in the (S] m , fi q ) plane, as in the LDGP 
case 0. 

The evolution of the Hubble parameter is given by 



E_ _ 3£(1 + zf [Q m + (1 + w q )fi q (l + z) iw « 
Ho~ 2(E+ TO") 



(14) 



This equation shows that E < for all z, as in the 
LDGP case, and that £ ^ as z ^ - 1 (equivalently, 
as a — > oo), for all w q > — 1. This behaviour is illus- 
trated in Fig. n Thus there is no super-acceleration, 
and no big rip singularity in QDGP models. The purely 
gravitational phantom behaviour does not produce the 
pathologies associated with phantom matter in GR. 

Since E is always negative, the deceleration parameter, 
q = —(1 + E/HqE 2 ), can be never less than —1, as seen 
in Fig. n Note that the asymptotic future value of q 
depends on w q : 



3(1 + w q ) — 1 as z 



(15) 



The transition redshift z a at which decelerated expansion 
turns into accelerated expansion, is defined by q — 0, i.e., 



E{z & ) = -H E{z & f . 
For the parameter values in Fig.^ z & ~ 0.7. 



(16) 
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redshift 

FIG. 2: The equivalent GR phantom potential Ves(z) for the 
parameter values in Fig. 
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FIG. 3: Phase plane trajectories for w q = —0.8. 



III. DYNAMICAL SYSTEM ANALYSIS 



The QDGP and LDGP models show an important dif- 
ference in the asymptotic behaviour of the Hubble rate: 



E 



0, 



-1, 

-1 . 



(17) 



In LDGP models, the de Sitter solution is a stable at- 
tractor. This arises because the vacuum energy does not 
redshift. By contrast, in QDGP models we have H — > 0, 
since the quintessence is redshifting away. 

The quantities p e ff and w e g describe a phantom scalar 
field c ff in GR which has the equivalent expansion his- 
tory to the QDGP model. We can find the self-interaction 
potential V e f[((j) c g) of the equivalent GR phantom field 
by using the well-known result that a scalar field and 
its potential may be determined for any given expansion 
history as follows 0: 



Vcs(z) 

<ftcff(z) 

V3 



dz 



(l + z)E 



3 _ d[E 2 - n m (i + zf 

2 dln(l + z) 3 

d[E 2 - n m (i + z) 3 ] 

dln(l + z) 3 



(18) 



1/2 



.(19) 



In principle, these equations lead to V e s — V c g(4> c g), but 
in practice the inversion cannot be performed analytically 
because of the integral in Eq. (|19fl . In Fig. [2] we show 
V e g(z). For the LDGP model, V e g approaches a constant 
nonzero value as a — > oo. By contrast, in the QDGP 
model the potential vanishes asymptotically. 



In order to write the cosmological equations in the form 
of an autonomous system, we define the following expan- 
sion normalized variables: 



" a 3 / 2 E ' V a 3(i+^)/2 E 
so that Eq. I|12H becomes 



E 



2v = x l 



y 



l 



(20) 



(21) 



This constraint means that the phase space is defined by 
x +y 2 > 1, since v > 0, and by x > 1, y > 1. Introducing 
the new time variable r = In a, and eliminating v and E, 
we obtain the autonomous system 

, _ 3a; [x 2 + (1 + 2w q )y 2 - l] 



2(x 2 



1) 



3y [2x 2 + (1 + w q )(y 2 - x 2 - 1)] 
2(x 2 +y 2 + 1) 



(22) 

// = — — (23) 

The analysis depends on whether w q > -1 or w q = -1. 

For w q > —1, the phase space is 2D. The fixed points 
and their stability are summarized in Tabled In this case 

TABLE I: Fixed points for > w q > -1. 
Point (x, y) Eigenvalues Character 



M (1,0) (-5=i,|) 



repellor 
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TABLE II: Fixed points for w q = -1. 



Point 



Eigenvalues Character 



M 



ds ^/n r jn A + v^ + ajav 



repellor 



attractor 



10 



20 



y 



30 



40 r dS 



FIG. 4: Phase line trajectories for w q = — 1, y/Q rc /^A = 21. 



there are two isolated fixed points located at finite values 
of x and Combining this fact with a simple examina- 
tion of the y' equation and the open topology of the phase 
space, we can draw fairly general conclusions about the 
properties of the fixed points of the system even before 
doing the linear analysis. Since y > and w q > — 1, we 
have y' > 0, so that the fixed points cannot be stable. 
In addition, the fixed points cannot be unstable spirals 
or centres, so they can only be unstable nodes (attrac- 
tors) or saddle points. The linear analysis confirms these 
findings. The fixed points M (early time) and Q (late 
time) are matter- and quintessence-dominated solutions 
respectively. The phase plane trajectories are illustrated 
in Fig. 

For the LDGP limit, w q = — 1, we have y oc v, leading 
to a ID phase space, and the dynamical equations reduce 
to 



3y(l-y 2 + 2 y /n r jn A y 



2(l + y/nr c /Q A y 



(24) 



There are also two isolated fixed points at finite values 
of y (see Table Inland Fig.0J. It can be seen that y' > 
as well, but since the phase space is only a line segment, 
and the fixed points are just its end points, it follows that 
one fixed point will necessarily be unstable (the matter- 
dominated point M), and the other will be forced to be 
stable (the de Sitter attractor dS). 



IV. CROSSING THE PHANTOM DIVIDE 

Perhaps the most distinctive feature of the QDGP 
model is the possibility of crossing the phantom divide 
w = — 1. This occurs with non-phantom quintessence, 
via the infrared gravitational modifications to GR. In 
the standard GR model, this can be realized only by 



introducing additional complicated features like multiple 
dark energy fields or interactions between dark energy 
and CDM. 

Equation J7J defines the effective equation of state, 
which leads to 



1 + w cS (z) = 

(i + w q )n q E(z)(i + z )3«»,+3 



.(25) 



/OT] + z)3» q +3 _ 2y^~E(2 

For the LDGP model, i.e., the io q = -1 limit of QDGP, 
w e g is always more negative than — 1 , provided that p e g > 
0. When p c g = 0, at some redshift z*, we have w e s = 
— oo. For z > z t , we have w g > — 1, with w e s — * +co as 
z — > z*. Thus w e g never passes through — 1. 

For w q > —1, the nonzero term (1 + u; q )O q £'(^)(l + 
z) 3u,q+3 in Eq. I|25|l allows for w e s to pass through — 1. 
This is illustrated in Fig. O For w q > —1, smooth cross- 
ing of the phantom divide occurs at a redshift z c that 
depends on the values of the free parameters f2 m , f2 q and 



By Eq. z c is given by 



(i + z c y w -E{ Zc ) 



(1 + W q )Q q 



(26) 



Note that for LDGP (w q = —1), the solution is z c = — 1, 
i.e., crossing occurs at a = oo. 

At a redshift z* > z c , we have p e ff(z*) = 0, so that the 
effective phantom GR picture of QDGP breaks down: 



4 j 1 ^""^ 



l/3(21+t«q) 



1 



(27) 



As in the LDGP case, we have w c s — > ^Fco as z 




FIG. 5: The effective dark energy equation of state w e ff(z) 
for the same parameter values as in Fig. 
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V. CONCLUSIONS 

The DGP(— ) class of braneworld models can lead to 
phantom-like behaviour of the effective dark energy, but 
without the need for any phantom matter, as pointed out 
by Sahni and Shtanov Q ■ In the simplest LDGP model, 
the universe contains only cold dark matter and a cos- 
mological constant. We generalized this by introducing 
a quintessence field, to define the QDGP model. QDGP 
has the same effective phantom behaviour as LDGP. The 
Hubble parameter is always decreasing, and there is no 
big rip singularity in the future. The key new feature of 
QDGP is crossing of the phantom divide, w = — 1, which 
is not possible in the LDGP case. 

The avoidance of any big rip is due to H, H — > as 
a — ► oo. This asymptotic behaviour reflects the fact that 
the phantom effects never dominate - unlike the case of 
a phantom scalar field in GR. The total equation of state 
parameter, defined by u> to t = Ptot/ptot = w e gp e $/{p + 



Peg), follows from Eqs. lfT2"jl and lj2"5)l. 

1 + Wtot(z) = 

n m (i + zf + (i + w q )n q (i + z)3(i+^) 

E(z)[y/n^ + E{z)] 
This shows that w e g(z) > — 1. 
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